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Abstract. We demonstrate that backward stochastic differential equations 
(BSDE) may be reformulated as ordinary functional differential equations on 
certain path spaces. In this framework, neither Ito's integrals nor martingale 
representation formulate are needed. This approach provides new tools for 
the study of BSDE, and is particularly useful for the study of BSDE with 
partial information. The approach allows us to study the following type of 
backward stochastic differential equations 

d 

dY^ = -f^{t, YuL{M)t)dt - J2 f^(t> Yt)dBl + dMl 

1=1 

with Yt = ^, on a general filtered probability space {Vt,F.Ft,P), where B 
is a c?-dimensional Brownian motion, L is a prescribed (non-linear) mapping 
which sends a square-integrable M to an adapted process L{M), and M, a 
correction term, is a square-integrable martingale to be determined. Under 
certain technical conditions, we prove that the system admits a unique solu- 
tion (Y, M). In general, the associated partial differential equations are not 
only non-linear, but also may be non-local and involve integral operators. 
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1 Introduction 



Stochastic differential equations (SDE) may be considered as dynamical systems 
perturbed by random signals which are often modelled by Brownian motion. 
The important class of stochastic differential equations considered in literature 
are Ito's type equations such as 

d 

dXi = }l{t,Xt)dt + ^ fl{t,Xt)dBl, (1.1) 
1=1 

where B ~ {B^, • • • , i?"*) is Brownian motion in on a completed probability 

space {Q,J-,P), fi = X]^=i fi'^ ^^"^ bounded, smooth vector fields in , 
j = 1, - ■ ■ ,d', where d, d' are two positive integers. K. Ito gave the meaning of 
solutions to (|1.1[) by developing a theory of stochastic integration against Brow- 
nian motion (called Ito's calculus), and obtained strong solutions by specifying 
an initial data at a starting time T. 

SDE (|l.ip has to be interpreted as an integral equation 

Jq Jo 

which can be solved forward (i.e. for t > 0). Ito's calculus requires that a 
solution X = (Xt) has to be adapted to Brownian motion B = (i?^, • • • , B"^), 
it is thus not necessarily possible to solve (jl.ip backward from a certain time T 
tot<T. 

There are interesting applications on the other hand to be able to solve 
(|l.ip backwards. Suppose m is a smooth solution to the Cauchy problem of the 
quasi-linear parabolic equation 

with the initial data u{x,0) — uq{x). Let T > and h(t,x) — u(T — t,x) for 
t e [0,T]. Then h solves the backward parabolic equation 

^A + ^^h + f{h,Wh)^0 on[0,T]xR'' 

and h{x,T) = uo{x). Let Yt — h{t,Bt) where B is Brownian motion in R'^. 
According to Ito's formula 

^T-yt=J^ (^-^ + ^A^h{s,B,)ds + MT-Mt (1.2) 

for i < T, where Mt — V/i(s, Bs).dBs is a martingale. Substituting + ^A) h 
by -fo{h,Vh) in pT^ to obtain 



Yr-Yt^-J^ f{Ys,\/h{s, Bs))ds + Mt - M*. 



(1.3) 
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According to Ito's martingale representation theorem, the density process Zt = 
\/h(t, Bt) of M with respect to Brownian motion is uniquely determined as the 
unique predictable process Zt such that 

Mt-Mo = V r ZldBl. 

3 = 1 ^0 

In terms of the pair (Y, Z) (|1.3p may be written as 

Yr-Yt^^ r f{Y,, Zs)ds + Y.f ^sdBi 

with the terminal data Yt — uo(-Bt), which is the integral form of the following 
backward stochastic differential equation 

dYt = -f{YuZt)dt + Zt.dBt, Yt^^ (1.4) 

introduced and studied in Pardoux-Peng |33j . 

In the past twenty years, there is a large number of articles devoted to the 
theory of BSDE and its applications in various research areas. Our references 
listed at the end of the paper is by no means complete, and the reader should 
refer to excellent surveys such as articles in the book jTH] edited by El Karoui 
and Mazliak, the recent paper by El Karoui et al [TJ and the references therein 
for a guide to the BSDE literature. 

To the knowledge of the present authors, it was Bismut ^ (see [7], [8]) who 
first formulated terminal problems for a class of stochastic differential equations 
in order to study stochastic optimal control problems by means of Pontrya- 
gin's maximum principal. His equations, called backward stochastic differential 
equations, have been extended and developed to a non-linear case in the seminal 
paper [33| by Pardoux and Peng. A lot of efforts have been made to generalize 
the class of BSDE considered in [35^. For example, Lepeltier and San Martin 
pS] relaxed the Lipschitz continuous condition on the driver and studied BSDEs 
with coefficients of linear growth. Yong [32] employed the continuity method 
to prove the existence of solution with arbitrary time horizon. In [9| Briand et 
al. considered L^-solutions for BSDE. It is also natural to consider BSDE cou- 
pled with a forward stochastic differential equation, called a forward-backward 
stochastic differential equations. Antonelli [1] first studied such FBSDE, his 
equation doesn't involve a density process Z in the driver. A definite account 
about FBSDE may be found in Ma et la [29], Hu and Peng [25], Peng and 
Wu [37], the recent book [30 and the literature therein. Most authors con- 
sider BSDE on a probability space with Brownian filtration, and there are a 
few papers dealing with BSDE with jumps or with reflecting boundary condi- 
tions. Tang and Li PT| have studied BSDE with random jumps, and Barles et al 
[4] have explored the connection between BSDE with random jumps and some 
parabolic integro-differential equations. Rong [39] proved the existence and 
uniqueness under non-Lipschitz continuous coefhcients for this class of BSDE. 



3 



Analogous to free-boundary PDE problems, El Karoui et al [TH] introduced an 
obstacle to BSDE such that the solution always stays above such obstacle. This 
so-called reflected BSDE is further developed to double reflected barriers by 
Cvitanic and Karatzas [TS], and Hamadene et al [53]. Furthermore, Barles et 
al |3] have considered BSDE on the probability space associated with Dirichlet 
processes. 

If the driver of BSDE is with quadratic growth of Z , the nature of equations 
is completely changed. This problem is firstly solved by Kobylanski [27j by using 
the monotonicity method adopted from the PDE theory. Her results have been 
substantially developed and generalized by Briand and Hu [TD] , [H] , where they 
extend to equations with convex drivers subject to unbounded terminal values. 
Most of the existing literature concentrates on solutions of BSDEs in a strong 
sense, i.e. the underlying filtered probability space is given. One of the first 
attempts to introduce weak solutions for BSDEs was presented in Buckdahn et 
al [12] , and Buckdahn and Engelbert [T^ further proved the uniqueness of their 
weak solutions, while the coefficients of their BSDEs do not evolve a density 
process Z. On the other hand, the notion of weak solution for FBSDEs was 
introduced by Antonelli and Ma [2] and further developed by Ma et al [31] by 
employing the martingale problem approach. 

The backward stochastic differential equations have found many connections 
with other research areas: stochastic control, PDE, mathematical finance and 
etc. To derive a maximum principle as necessary conditions for optimal control 
problems, one can observe that the adjoint equations to the optimal control 
problems satisfy certain backward equations. For stochastic control problems, 
the corresponding adjoint equations are stochastic rather than deterministic. 
Indeed Peng [35] established a general stochastic maximum principle by consid- 
ering both first order and second order adjoint equations, and, on the other hand, 
Kohlmann and Zhou [26] interpreted BSDE as equivalent to stochastic control 
problems. Peng |36| derived a probabilistic representation (a Feynman-Kac 
representation) for solutions of some quasi-linear PDEs, which was extended to 
other cases by Ma et al [5S]. The later has been summarized as a four-step 
scheme of solving forward-backward stochastic differential equations (FBSDE) , 
see [30] by Ma and Yong for a detail. In [TB] DuflSe and Epstein discovered a 
class of non-linear BSDE in their study of recursive utility in economics. Later 
El Karoui et al [21, applied BSDE to option pricing problems and provided a 
general framework for the application of BSDE in finance. In order to deal with 
utility maximization problems in incomplete markets. Rouge and El Karoui j40j 
introduced a class of BSDE with quadratic growth. Hu et al [24j further studied 
this class of BSDE in a more general setting. 

In this article, we put forward a simple approach to deal with the kind of 
BSDE such as (|1.4p which does not depend on any martingale representation, 
thus allows to study a wide class of backward stochastic dynamics. Our main 
idea and contribution in this article is to establish an ordinary functional differ- 
ential equation which is equivalent to (|1.4p . which allows us to obtain alternative 
representations for solutions of BSDE and to consider a new interesting class of 
stochastic dynamical systems. 
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Consider the following example of backward stochastic differential equations 



dYt = -/(t, Yt, Zt)dt + ZldBl Yt = 



(1.5) 



where B — {Bt)t>o is Brownian motion in R'^ , ^ G L'^i^, Ft, P), and {Ft)t>o is 
the Brownian filtration. The differential equation has to be interpreted as the 
integral equation 



Jt ._1 Jt 



(1.6) 



By applying the Picard iteration to (Y, Z) , one shows that if / is globally Lip- 
schitz continuous, then there is a unique pair (Y, Z) which satisfies (11.61) for 
all t < T. This method relies on the martingale representation for Brownian 
motion, and thus restricts the class of BSDE. 

Our main idea is based on the following simple observation. Suppose that 
Y = (Ff )tg[^ is a solution of (II. 6p back to time t < T, then Y must be a special 
semimartingale whose variation part is continuous. Let Yt — Alt — Vt be the 
Doob-Meyer's decomposition into its martingale part M and its finite variation 
part —V. The decomposition over [t,T] is unique up to a random variable 
measurable with respect to J>. Let us assume that the local martingale part 
M is indeed a martingale up to T. Then, since the terminal value Ij- = ^ is 
given, ^ = Mt - Vt, so that Mt = E{C + VT\Ft) and Yt = + Frl^t) - V 
for t G [t^T]. The integral equation (II. 6p thus can be written as 



i-T d „T 

^-Mt + Vt^- f{s,Ys,Zs)ds + J2 Z. 

Jt Jt 



dBl 



for every t G ['T',?^]. 
obtains 

m\^t) 



Taking expectations both sides conditional on Ft one 



Mt + Vt = 



fis,Y,,Z,)ds 



Ft 



fis,Y.s,Zs)ds. 

By identifying the martingale parts and variational parts, we must have 



Vt-Vr 



f{s,Ys,Zs)ds 



where Y and Z are considered as functionals of V, namely 

Yt = E{^ + VT\Ft) - Vt, Mt = EiC + VT\Ft) 



(1.7) 



(1.8) 



5 



and Z is determined uniquely by the martingale representation through 

d „T 

Mt - Afr = ^ / Z\dB\. 
1=1 "'^ 

Hence Y and Z are written as Y{y) and Ziy") respectively, if we wish to em- 
phasize the fact that Y and Z are defined entirely through V . Observe that 
(|1.7p is clearly the integral form of the functional differential equation 

which can be solved by Picard iteration applying to Y alone, rather than the 
pair {Y,Z). 

The approach may be made independent of the use of a martingale repre- 
sentation theorem, provide that one is willing to replace the density process Z 
by a functional of Y , thus free us from the requirement of Brownian filtration. 
This kind of generalizations of BSDE theory comes a little bit surprising and 
even over rewarded, which is though not the only point we would like to empha- 
size. More precisely, we may consider the correction martingale part appearing 
in (|1.5|) as part of the solution rather than its density process Z. That is, by 
setting Mt — Mt = Ylt=i It ^I'^^l^ ^^'-^ regarding Z as a function of M, so 
denoted by L{M), then p. 51) can be reformulated as 

dYt^~f{t,YuL{M)t)dt + dMt, YT = i, (1.9) 

which is in turn equivalent to the functional integral equation 



Yt-YT = 




f{s,Y{Y)s,L{M{Y))s)ds (1.10) 



where Y{Y) and M{Y) are given by l|1.8p . For (jl.lOp . there is no need to insist 
that L sends a martingale AI to its density process (if there is any), though the 
density process mapping L remains the most interesting case. 

The approach might be applied to a more general setting of solving dynam- 
ical systems backward under other constraints, not necessarily the adaptedness 
to a filtration, even a probability setting is not necessary. One possible example 
can be the following. One may study the functional differential equation (|1.7p . 
where Y : Y Y{Y) and M : Y ^ MiV) are defined in terms of some kind 
of "projections" instead of conditional expectations. We however in this paper 
make no attempt for such an extension. 

To our knowledge, most of BSDE which currently exist in the literature may 
be studied in the framework of ordinary functional differential equations. Since 
our approach does not rely on the martingale representation theorem, we are 
able to study a class of BSDE on an arbitrary filtered probability space. We 
however would like to point out that this paper is not so much about gener- 
alizing the theory of BSDE to a general filtered probability space, our main 
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contribution is the equivalence of BSDE and a class of ordinary functional inte- 
gral equations. We allow a sufficient wide class of functionals L{M) which even 
in the classical setting extends the associated PDE to some no-local integro- 
differential equations. 

If {J-t)t>o is Brownian filtration, any martingale M = {Mt)t>o has an Ito's 
integral representation Mt — Mq = J2'j=i lo ^sd-W^ which determines the den- 
sity Z — {Z^, • • • , Z"^). Consider the functional over martingales 

where fj,(t,ds) is a transition kernel (not random for simplicity), and consider 
the corresponding BSDE 

dYt' = -f^iT - t, YuL{M)t)dt -f dMl Yt = 

Our approach applies to demonstrate the existence and uniqueness for this kind 
of BSDE, whose associated PDE is a system of integro-differential equations: 

where the non-linear operator H involves space-time integration, and indeed 

H{u){t,x)= / —j -i^~^)dzti{t,ds). 

Jt Jr" (27r(s-t))''/^ 

If fi{t, ds) — St{ds) then we recover the case considered in the current literature. 
By choosing different functional L{M) we may obtain even more general integro- 
differential equations. This kind of integro-differential equations often appear 
in the study of particle limiting models for PDE, one class of equations which 
has the similar nature are already in the literature, for example in Majda [35]. 

In this paper we constrain ourself to the study of the following type of 
backward stochastic differential equations 

d 

dV/ = -f^{t, Yt,L{M)t)dt - ^*)^^? + (1-11) 

1=1 

subject to Yt = on a filtered probability space {n,J^,J^t,P), where i? is a 
(i-dimensional Brownian motion as given, j = 1, • • • , d', L is a given (non-linear) 
functional on square- integrable martingales, while {J^t)t>o is not necessary to 
be Brownian filtration. A solution to (|l.Iip is a pair {Y,M), where Y = (Y^) 
are semimartingales and M — [A'P] are square-integrable martingales, which 
satisfies the corresponding integral equations: 

Y^ - Yl ^~ r f^it, Yu L{M)t)dt -Y. r flit, Yt)dBl + - Mi. (1.12) 

-'t ^=l •'t 
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The term L{M) appearing in the drift term /o on the right-hand side of 
suggests that L is a mapping which sends a vector of square-integrable martin- 
gales M = (M^) to a progressively measurable process L{M). The backward 
stochastic equation is thus described by the driver /o, the diffusion coef- 

ficients fi together with the prescribed mapping L. 

Finally, let us point out that similar ideas have been known in the PDE 
theory. Recall that, for any reasonable function u, u has the following decom- 
position: 

u = H{u) + G{u) 

where H{u) is a harmonic function determined by a boundary integral against 
a Green function, and G{u) is a potential. Thus the boundary condition (which 
corresponds to our case the terminal value) determines the harmonic function 
part H{u). The regularity theory for non-linear PDE such as Am = /(m, Vu) 
may be developed via the previous decomposition, by studying the Newtonian 
potential G{u), (Gilbarg and Trudinger [12]). In this way, backward stochastic 
dynamics, as a class of Markov processes, can be regarded as a generic exten- 
sion of some non-linear PDE problems of finite dimension to infinite dimensional 
problems in path spaces. On the other hand, some non-linear PDE can be con- 
sidered as a pathwise version of backward stochastic dynamics. We will explore 
these ideas further in coming papers. 

The paper is organized as following. In Section 2 we present some elementary 
facts and basic assumptions. The main results of the existence of local and 
global solutions, and the uniqueness of the backward stochastic dynamics are 
presented and proved in Sections 3 and 4. 



2 Preliminaries 

Let (f2, J^, Tt-, P) (where t e [0, oo)) be a filtered probability space which satisfies 
the usual conditions: {il,J-,P) is a completed probability space, (J-'t)(>Q is a 
right-continuous filtration, each J^t contains all events in J- with probability zero, 
and J-' — (7{J-'t : t > 0}. Under these technical assumptions, any martingale on 
{Q,J^,J^t,P) lias a modification whose sample paths are right continuous with 
left-hand limits. Henceforth, by a martingale we always mean a martingale 
which is right continuous with left-hand limits. 

Let < r < T be any but fixed numbers. [t,T] serves as the region of 
the time parameter, although we are working with a fixed filtered probabil- 
ity space {fl,T,Tt,P)- Let C ^[t,T] denote the space of all continuous, 

adapted processes {Vt)te[T.T] valued in J?'* such that max^ supjgj^ |V/ 1 belongs 
to L^{^1, J-T, P), equipped with the norm 



C[r,T] 



J2e sup \Vt 
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C \\t,T] ;R'^j is a Banach space under || • ||c[r,T]- M'^ {[t , T]; R^) denotes the 

space of fi''- valued square-integrable martingales on (VL,!F ^TtjP) from time r 
up to time T (which, of course, can be uniquely extended to a martingale in 
M'^{p,T],R'^)), together with the norm ||M||c[r,T]- We also need the direct 

sum space of M"^ {[t ,T]\ R"^) and C {[t,T];R'^^, denoted by S{[t,T]\R^). If 
Y e 5([t, T]; i?''), its decomposition into an element in A^^([r, T]; ii'') and the 
other in C ^[r,T] ; .R'^^ may not be unique, and there are various norms one 
can define on S{[t ,T]- R"^). For our purpose, we choose the norm ||y||c[T.T]j 
although 5([t, T]-R'^) is not complete under || -Hcir.T]- Finally ^^([r, T]]R'^''"^) 
denotes the space of all predictable processes Z = {Zf'^)t(z[r,T] on {fl, J^, J^t, P) 
with running time [tjT], endowed with the usual i^-norm 

d' d „T 

j=i -^^ 

If y is a semimartingale on {il,JF,J^t,P) over time interval [t,T] with its 
Doob-Meyer's decomposition Yt = Mj — Vt, such that M is an J't-martingale 
during [r,T], y is a continuous, adapted process with finite variation on [T,r], 
and Vt, Yt are integrable, then Mt = E{YT + VT\Tt) and Yt = E{YT + VT\Tt)- 
Vt for t € [t,T]. Since we are interested in terminal value problems, in which 
Yt — £, are given, therefore, for given ^ = (^*) where C & L^i^i^T^P), we 
consider two functionals on C {[t, T] ; R'^^ : V Y{V) and V M (V) defined 

by 

Y{V)t^E{^ + VT\J't)-Vt forte[T,r] (2.1) 

and 

M{V)t= E{^ + VT\Tt) for fe[T,T] (2.2) 

for any V G C (^[t , T] ; R'^^ . If we wish to indicate the dependence on the 

terminal value ^ as well, then we will use Y{£^, V) and M(^, V) in places of 
Y{V) and M{V) respectively. 

Note that (Y{V)t)t£iT,T] does not depend on the initial value Vr, an im- 
portant fact we will use in our construction of global solutions for the terminal 
value problem ()l.lip . 

We consider the following type of backward stochastic differential equations: 

d 

dY^ = -fi{t, Yt,L{M)t)dt - Y fiit, Yt)dBl + dM/, Y^ = ^ (2.3) 

on the filtered probability space (fi, J^, P) (j = I,-- - ,d'), where i? is a 
d-dimensional Brownian motion on ,J-t,P) as given, T > is the ter- 

minal time, e {Q, , Ft , P) (for j = I,-- - ,d') are terminal values, // 
(z = 0, • • • , d and j = 1, - ■ ■ ,d') are locally bounded and Borel measurable, 



\n' 



\ 
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and L is a prescribed mapping on M'^{[t,T]; R!^' ) valued in H'^{[t,T]; R!^' '''^) 
or in C ([T,T];R'^'y 

A solution of (|2 3p backward to time r is a pair of adapted processes {Yt, AIt)t^[T,T] 
where = {Ml)tfz[r,T] are square-integrable martingales and 1^ = {Yf)i(z[r,T] 
are special semimartingales with continuous variation parts, which satisfies the 
integral equations 

r/ - = r fi {s, L{M),)ds + J2 r fi Y,)dBl + Mi - (2.4) 
Jt Jt 

forte[T,T],j = !,■■■ ,d'. 

As we have seen in Introduction, by writing Yt = Mt — Vj, a solution {Y, M) 
to p.4p is equivalent to a solution V of the functional integral equation 

Vt-V,= f Ms,Y{V)s,L{M{V)),)ds + J2 f Ms,Y{V),)dBl (2.5) 

where M(y)t = + FtI^*) and Y{V)t = M(y)t - Vt for t e [t,T]. It is the 
integral equation (|2.5p we are going to study. 

The following standard assumptions are always imposed on our backward 
SDE (|2.3p . Additional conditions on L will be introduced later on to ensure 
local and global existence. 

1) /o — {fi)j<d' are Lipschitz continuous on [0,oo) x R"^ x JJ™, and fi = 
{fi)j<d' (i = 1, ■ ■ • ,d) Lipschitz continuous on [0, oo) x : there is a constant 
C2 such that 

|/o(t,j/,^)| <C2(l + t+|y| + |z|), 
|/o(t, y, z) - h{t, v\ z')\ < C2i\y - y'\ + \z - z'\), 
m,y)\<C2il + t+\y\), 

and 

\Mt,y)^hit,y')\<C2\y-y'\ 
for t > 0, aU y, y' £ R'^' and z, z' e R"" . 

2) The terminal value £, = {C)i=i,- ,d' , C e L'^{n,TT,P). 

3 Local solutions and uniqueness 

In this section, we prove two results: the uniqueness and the existence of a local 
solution to ()2.3p under the assumption that L is Lipschitz continuous: 



3) L:MH[r,T];R''')^HH[T,T];R"') (resp. C{[t,T]; R"')): 
\\L{M) - L{M)\\n2 <Ci\\M-M\\c 
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(resp. 

\\L{M) - L{M)\\c < Ci\\M - M\\c) 

for any M, M e M^{[t,T]; R''-'), where ||M||c means ||M|lc([r T]-R<i') "^tc. for 
simplicity. 

By a local solution we mean a solution from T back to r, where T — r is 
smaller than a certain constant depending on L and // . 

In order to prove the uniqueness, we need to consider BSDE in a more 
general form than (|2.3p . More precisely, we are given another Brownian motion 
W = (T^i,--- on (f], J"t, J",P) and gk : R+ y~ which are 

Lipschitz continuous: 

\9k{t,y)\<C2{l + t+\y\) 

and 

\gk{t,y) - 9k{t,y')\ < Caly-y'l 
for &\\t>Q,y,y' eR'^' ,k^\,--- , m' . Define 

L:M^{[t,T];R''')xS([t,T];R''')^H^{[t,T];R"') (resp. C([r,T];fi™)) 

by 

L (M, r) = L I^A/ - E / 9k{s, Y,)dW^'^ . (3.1) 

Consider the following mapping L defined on Co{[t,T]]R'^'), the space of all 
processes in C([r, T]; R'^ ) with initial data Vr — 0, by 

L(v)t = ^ /o(s,r(y)„i(Af(y),r(y)),)ds 

+ f Ms,Y{V)s)dBi (3.2) 

where Af(y)t = + W| J"t) and y(V)t = M{V)t - Vt for t £ [r, T], so that 
Y{V)t = As we have seen that, the functional integral equation: V = L(T^), 
is equivalent to the following BSDE 

d 

dF/ = -f^{t, YuL{M, Y)t)dt - flit, Yt)dBl + dM^ Yt = ^ (3.3) 



Lemma 3.1 L defined by i3.1\) is Lipschitz continuous: 
\\L{M,Y)-L{M,Y)\\^.^,^T] 
< Ci||M-M||cKT] + ^^^^(T-r)||y-f||c[.,T] (3.4) 
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\L{M,Y)-L{M,Y)\\c[r,T] 



< Ci\\M - M\\c[r,T]+2m'CiC2VT^\\Y -Y\\c[r,T] 
for any M,M e M'^{[t,T]; R'^') and Y,Y € C{[t,T]; R^'). 

Proof. Let us omit the subscript [r, T] for simplicity. Then 
\\L{M,Y)-L{M,Y)\\n^ 



(3.5) 



< Ci\\M-M\\c + Ci 



fe=i 



Y,) - gk{s,Ys))dW^ 



= C^\\M - M\\c + Cif^J E £ j\gkis,Y, 



) - gk{s,Ys))dW^^ 



dt 



Ci\\M-M\\e + C,J2\^ f f 
fe=i V Jt Jt 



{gk(s,Ys) - gkis,Ys)) 



dsdt 



< Ci\\M - M\\c + m'CiC2\ E 



T rt 



dsdt 



< Ci||M-M||c + ^^^^(T-r)||F-y||c. 

The proof of the second inequaUty is similar. ■ 

The following is our basic local existence theorem. 

Theorem 3.2 Under the assumptions on L, /j and described above. Let 

1 



n = 



CI 



4Ci + 6 + 2Vdj + 3\/2m'CiC2 



A 1 



(3.6) 



which depends on the Lipschitz constants Ci,C2 and the dimensions, but is 
independent of the terminal data ^. Suppose that T — t < h, then L admits a 
unique fixed point on Co([t, T]; i?'' ). 

Proof. The proof is a standard use of the fixed point theorem applying to 
L. To this end, we need to show that L is a contraction on Co{[t,T]; R ) as 
long as T — T < h. This can be done by devising a priori estimates for L. Let 
us prove the case that L : M'^ {[t , T]; R'^' ) U'^ {[t ,T]; R"') is Lipschitz, the 
other case can be treated similarly. To simplify our notations, let 5 = T — Zi be 
the life duration. Since 



\\UV)\\c < ^y^/ \fo{s,Y{V)s,L{M{V),Y{V))s)\Hs 
+2 f \fi{s,Y{V),)\Hs. 
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Since /o and fi are Lipschitz continuous, so that 



\HV)\\c < 2C2 (v^ + Vd) (l + s)2ds 



+2C2 Vd) E\Y{V)s\^ds 

+2C2V6\\L{M{V),Y{Vmn2. 
Together with the elementary estimates 

\\Y{V)\\c<2./E\^ + 3\\V\\c 

and 



(3.7) 



||M(l/)||c<2^^ + 2||y||c 



one deduces that 



+ 2 \V2m'CiCid + 2C2V6 + 2C2\/rf + 2C2C1I V6^/E\^ 
+ [3V2to'CiC|(5 + 6C2\/^ + 4C2C1 +6C2\/dl \/^||F||c. (3.8) 
Similarly, for V,V G C[t,T] such that Vr = Vt = one has 



\\Uv)-uv)\\c < 



^E\^J \fo{s,Ys,L{M,Y)s) - h{s,Ys,L{M,Y)s)\ds 



\ E sup 
\ te[T,T\ 



E f [/.(5>n)-/i(s,n) 



dB\ 



where = E{^ + FriJ^t), = E{^ + FrlJ^t), = - Vt and = Mt - 14. 
Since fi are Lipschitz continuous, so that 



|/o(s, n, L(M, Y)s) - Ms, Ys, L{M, Y),)\ds 



< a 



E 



Ys-YA + 



■.{M,Y)s-L{M,Y)J^ds 



Ys-YJ + 



:{M,Y)s-L{M,Y)s\yds 



< C2S\\Y -Y\\c + C2VS\\L{M,Y) - L{M, y)||„2 



< C2 



1 + V5- 



m'CiC2 

71~ 



6\\Y-Y\\c + C2CiVS\\M-M\\c 
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where the last inequahty follows from p.Sp . The Ito's integration term can be 
treated similarly. Applying Doob's inequality, one has 



< 2 




dBl 



dBl 



< 2C2\/dv^||F-f||c. 
Putting these estimates together we obtain 



1 + %/^. 



m'CiCh 
V2 



m-Y\\c 



||L(^)-L(y)||c < C2 

+C2 (Ci + 2\/d) ^f5\\M ~ M\ 
On the other hand it is easy to see that 

\\M-M\\c = Je sup E{VT~VT\:FtY 

te[T,T] 



c- 



(3.9) 



< m~v\\c 



and 



\\Y - Y\\c < 3\\V - V\\c . 
Inserting these estimates into p.9p we finally obtain 

\\Uv)^Uv)\\c<C2 



2Ci+6Vd + 3V^+3^^^^^ 



VS\\V-V\\c. (3.10) 



Since 6 < li, the constant in front of the norm on the right-hand side is less 
than i, so that 

\\Uv)~uv)\\c<l\\v-v\\c. 

Therefore L is a contraction on Co([t, T]; R"^) as long as T — t < li, so there is 
a unique fixed point in Co[t,T]. This completes the proof. ■ 

We are now in a position to show the local existence and uniqueness of 
solutions to BSDE (lO). 



Theorem 3.3 Let L, /j be Lipschitz continuous with Lipschitz constants Ci, 
C2, and 

1 



k = 



4Ci + 6 (1 + 2Vd') + 3V2d'CiC2 



A 1 
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which is independent of the terminal data ^ G L'^(n,TT, P)- Suppose that 
T -T < h, and L{M) = L{M - Mr) for any M € {[t , T]; R'^' ) . Then 
there is a pair (Y,M), where Y — {Yt)i^[r,T] « special semimartingale, M = 
iMt)t£[T,T] is o, square-integrable martingale, which solves the backward stochas- 
tic differential equation \2. 3\] to time r. Moreover, such a pair of solution is 
unique in the sense that if (Y, M) and {Y , M) are two pairs of solutions, then 

Y ^Y and M - NU = M - Mr on [r, T] . 

Proof. By Theorem 13.21 (applying to the case that all — 0), there is a 
unique V G Co[t,T] such that 

fo{s,Ys,L{M)s)ds + Y, / f^{s,Ys)dBl, Vt e [t,T] 
1=1 "'^ 

where Mt ^ E {£, + VxlJ't) and = Mt -Vt- It is clear that Yt = i and 

Yt-^^ fo{s,Ys,L{M),)ds + Y, f^{s,Y,)dBl + Mt-MT (3.11) 
Jt ^^^Jt 

for all t G [t,T]. Therefore {Y,M) solves the backward stochastic differential 
equations (12.31) . 

Suppose that {Y,M) and {Y,M) are two solutions satisfying p. lip , where 

Y and Y are two special semimartingales. Let 



Zt = Mt+Y, f^{s,Ys)dBl 

,-1 J-r 



Then 



^t-^ = l fois,Ys,L{Z,Y),)ds + Zt-ZT (3.12) 
for t G [t,T], where 



L{Z,Y)^L\^Z-Y,J^ f,{s, Y,)dBl 
It follows that 

Yt=E[^ + AT\Ft]-At 



where 



At = fo{s,Ys,L{Z,Y),)ds, yt G [r,T]. 
Hence Yt = Y{A)t and the integral equation (I3.12[) becomes 

Yt^Ar-Zr + ^-At + Zt 
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Since Ar — so that 

Yr^Ar-Zr + ^ + Zr 
and thus we may rewrite the previous identity as 

Yt=Yr + [Zt - Zr) - At 

By the uniqueness of the decompositions for special semimartingales we must 
have 

Yr + {Zt -Zr)^E[i + AT\Ft] - M{A)t. 

Since L{M) = L{M - Mr) for any M e {[t ,T]\ R^' ) , so that L{Z,Y) = 
L{M{A),Y). Hence 

At = J' fo{s, Y{A)s,L{M{A), Y{A)),)ds. 

The same argument apphes to (Y, M), so that we also have 

At = J' fo{s, L{M{A), Y{A))s)ds. 

By Theorem [321 A = A, which yields that Y = Y. It follows then 

Zt - Zr = Zt - Zr yt e [t, T] 

thus M — Mr — M — Mr which completes the proof. ■ 

One of course wonders whether can the global existence be established, by 
means of weighted norms for example as in the BSDE literature. The present 
authors were unable to achieve better results than the local existence even with 
different choices of norms or spaces to which we apply the fixed point theorem. 
In fact, under the Lipschitz condition only on the mapping L, the local existence 
is the best we can hope, this is because L{M)t may depend on the whole path 
from r to T, and therefore the corresponding stochastic functional differential 
equation 

d 

dVt = h{t,Y{V)t,L{M{V))t)ds + Y,U{t,Y{V)t)dBl , K = 

1=1 

is neither local nor Markovian. This can be best demonstrated by its associated 
differential and integral equation. For example, it is not difficult to show that 

Lc{M)t = ^E{{M-,M-)t - {M-,M-)t\Ft) 

for t G [t, T] is Lipschitz continuous, where M'^ is its continuous martingale part 
such that Mq = 0, and therefore we have 
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Corollary 3.4 Suppose T <l2. Then there is a unique special semimartingale 
Y — {Yt)te[a,T] such that It = ^ o,nd 

fois,Y,,L,{M))ds + Mt-MT. (3.13) 
Moreover M is unique up to a random variable measurable with respect to Tq. 



Let us apply Corollarv l3.4l to the case that {J^t)t>o is the Brownian filtration 
of Brownian motion B = (B^ , ■ ■ ■ , B'^). Then, by Ito's martingale representation 
theorem, 



Lc{M)t 



\ 



T d 

Y,E{\ZmFt)ds 

i=i 



where are predictable processes such that 



mt - Mr f ^i^^i- 



Suppose u is a bounded, smooth function which solves the backward parabolic 
non-linear equation 

+ i + /o(i, u, K{u)) = on [r, T] x R'^ (3.14) 
with m(T, •) = (/?, where 



K{u){t,x) = ^ P,_t|Vw|2(s,x)ds 

where {Pt)t>o is the heat semi-group in R^, i.e. Pt = e^*^. In particular, the 
differential and integral equation p.l4p is not local, and is a non-linear equation 
involving space-time integration and partial derivatives. 

Applying Ito's formula to the process Yt = u(t, Bt) one has 

YT~Yt = ^ (^^^ + ^A^u{s,B,)ds + MT-Mt 

fo{s, Y,,K{u){s, B,))ds + Mt - Mt 

where Mt — /J Vu(s, Bs).dBs is a square-integrable martingale, and one recog- 
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nizes that 



Lc{M)t 



^E{{M,M)T-{M,M)t\Tt) 



\ 



\Vums,B,)ds\Tt 



P,_t|VM|2(s,Bt)ds 
= Kiu)it,Bt). 

Therefore {Y,M) is the unique solution to p.l3p . and we have a probabiUty 
representation 

uit,x) = E{Yt\Bt^x}. 

Since the non-hnear equation p.l4p depends on the "future" of the solution 
from time T, it is not always possible that a solution exists back to any time 
T. In turn, we thus can not expect that the general BSDE (|2.3p have a solution 
global in time without further restrictions on L. 



4 Global solutions 

In the previous section, under only the Lipschitz conditions on L we are able to 
construct a solution to the backward stochastic differential equation (|2.3|) back 
to time r such that T — t < 12- 

In this section we construct the unique global solution to (|2.3p if L satisfies 
further regularity conditions. 

We assume that the mapping L : M^{[0,T]; R!^') n^{[0,T]; R"") (resp. 
C([0, T]; ii™)) satisfies three technical conditions a), b) and c) below: the local- 
in-time property, the differential property, and the Lipschitz condition. The last 
one is standard, but the first two properties are motivated by the example of 
density processes in Ito's martingale representations. 

For any [T2,Ti] C [0,r], define the restriction 

L[T,,T,]-M\[T2,Ti];R''')^H\[T2,Ti];R"') (resp. C([r2, Ti]; fi"*')) 

by L[T2,Ti](^)t = L{N)t for any N e M^{[T2,Ti]; R'^' ) and t e [T2,Ti], where 
N e M^{[0,T]]R'^') defined by Nt = E{NT^\Tt) for t < Ti and Nt = Nt, for 
t > Ti. 

a) (Local-in-time property) For every pair of non-negative rational numbers 
T2 < Ti < T, and for any M G M'^{[Q,T];R'^' ), L{M) = L[t^^t,]{M) on 
(T2,ri), where M = (A/t)te[T2,Ti] is restriction of M on [T2,Ti]. The local-in- 
time property requires that L{M)t is locally defined, i.e. L{M)t depends only 
on {Ms)s^[t^t+E) for whatever how small the e > 0. 
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b) {Differential property) For every pair of non-negative rational numbers 
Ti < Ta < T, and M € M^{[T2,Ti]; R'^'), one has L[T2,Ti](M - MyJ = 
^{T2,Ti]{M) on (T2,Ti). The differential property requires that ^[j'j^Ti] 
depends only on the increments {Mg — Mt^ '■ s>t} for f e [T2,Ti]. 

c) {Lipschitz continuity) L : M'^{[0, T];R'^') "H^qq, T]; il™) (rcsp. C([0, T];R"')) 
is bounded and Lipschitz continuous: there is a constant Ci such that 

||L(M)||^.^^_^^^ <Ci||M||c[T„T,] (4.1) 

and 

||L(M)-L(M)||„.^^_^^^ <Ci||M-M||c[T„T,] (4.2) 

(resp. 

\\L{M)\\c[T„T,] < Ci\\M\\c[T„T,] (4.3) 

and 

\\L{M) - L(M)i|c[T„T,] < CilIM - M||c[T„T,]) (4.4) 

for any M,M G M'^{[0,T];R'^') and for any rationales Ti and T2 such that 
< T2 <Ti < T. That is to say L[X2.Ti] are Lipschitz continuous with Lipschitz 

constant independent of [T2,Ti] C [0,T]. 

The first example below provides the most interesting examples of L in ap- 
plications, which are however variations of the classical example considered in 
the literature. 

Example 1. Suppose that (J-t)f>o is the Brownian filtration generated by 
a d + A/'-dimensional Brownian motion B = {B^,- ■ ■ , B'^, , ■ ■ ■ , W^) on a 
probability space {n,T,P). If M e M'^{[0,T]; R"^'), then, according to Ito's 
martingale representation theorem, M is continuous, and there are unique pre- 
dictable processes (Zj'')tg[o,T] such that 

Mi -Mi = Y, Zi''dBi + ^ / Zi'^+'^dW^, j = l,- ■■,<!' (4.5) 

for ah t e [0,T]. Assign M G M^{[0,T];R'^') with L{M) = {Z^'')j<d',i<d- For 
< T2 < Ti < T, the restriction of M on [T2, Ti], denoted again by M, belongs 
toM.'^{\T2,Ti\;R'^ ). By the uniqueness of Ito's representation we can see that 
L satisfies the local-in-time and differential properties. It is also easy to show 
that L : M'^{[Q,T];R'^') 'H'^{[0,T];R'^''"^) satisfies the Lipschitz condition. 

Another class of interesting examples of L is presented in the following ex- 
ample. 

Example 2. Let {Q,J^,J^t,P) be a filtered probability space which satisfies 
the technical conditions described at the beginning of section 2, but not neces- 
sary to be a Brownian filtration. Let B = {Bt)t>o be a Brownian motion in i?"* 
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adapted to (J-t)t>o, therefore (Jt)t>o is in general bigger than the Brownian 
filtration generated by B. Let A4b denote the closed stable sub-space of A^2 
determined by B, that is 



T 



TUs = l^y" ZidBl : are predictable and E J 
Then any martingale M has a unique decomposition 

Mt - Mo = V / zidBi + m; 

where M' e M'^{[0,T]; R) orthogonal to Mb- Then L(M) = [Z^) satisfies the 
local-in-time and differential properties, as well as the Lipschitz condition. 

In the following theorems we retain the basic assumptions on the coefficients 
/j and the terminal values ^\ 

Theorem 4.1 Assume that L satisfies the conditions a), b) and c) listed above. 
Then there exists a pair of processes (y,M), where Y = (i^t)te[o.T] special 
semimartingale, and M = (-^t)te[o,T] a square integrable martingale, which 
solves the backward equation 

d 

dYt = -fo{t, Yt,L{M)t)dt - J2 Mt^ yt)dBl + dMt, Yt = (4-6) 

1=1 

The solution Y is unique, its martingale correction term M is unique up to a 
random variable measurable with respect to J-q. 

The remainder of this section is devoted to the proof of Theorem 14. II Recall 
that 

k = z -2 A 1 



4Ci + 6 ( 1 + 2Vd) + 3V2d'CiC2 



which is positive and independent of 

By Theorem 13.21 if the terminal time T < I2, the non- linear mapping L on 
Co([0,r];ii'^') admits a unique fixed point, where 



[ fois,YiV),,LiMiV))s)ds + J2 [ f^{s,YiV),)dBl 
Jo ,_i Jo 



UV) 



Next we consider the case T > k- In this case we divide the interval [0,r] 
into subintervals with length not exceeding Z2. More precisely, let 

T = To > Ti > • • • > Tfe = 
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so that < Ti-i — Ti < I2 where Ti are rationales except To = T. 

Begin with the top interval [Ti , Tq] , together with the terminal value Yto = ^ 
and the filtration starting from J-t^ . Applying Lemma [3.2l to the interval [Ti, Tq] 
and Li, where 

{UV)t = f Ms,Y,{V)s,L[T,,To]{Mi{V))s)ds 
+ E f U{s.Y^{V)s)dBl 

where 

M^{V)t = E{i + FtoI , Y^{V)t = h'h{V)t - Vt 

for any V £ C{[Ti,Tq]; R"^' ) and t € [Ti,ro]. Then, there exists a unique 
y(l) e Co{[Ti,Tq];R'^') such that LiV^(l) = V{1), 

Repeat the same argument to each interval [Tj,Tj_i] (for 2 < j < k) with 
the terminal value Yj-i{V{j — 1))tj_i , the filtration starting from , and the 

non-linear mapping Lj defined on Co([Tj, Tj_i]; i?'^ ) by 

mV)t = f fo{s,Y,iVUL[T„T,..]iMiV,))s)ds 

where V € C([T, , T,-_i]; R^') and 

M,{V)t - E{Y,^,{V{j - 1))t,_, + V^T,_J^t), 

fort e [r„T,_i]. 

Therefore, for 1 < j < fc, there exists a unique e C([Tj , T,_i]; il'*') 

such that 

= r/o(s,r(j%,£[T,,T,_,](M(j))s)ds 

for t e where Y{0)t„ = L Y{j ~ 1)t,_, - Y{j)t^_, for 2 < j < k, 

and 

M(j)t = - + y(j)T,_jJ-t), 
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for t e [Tj,T,_i]. 

Since y(i - 1)t,_i = Y{j)t^_, for 2 < j < fc, F = (yt)te[o,T] given by 

Yt = Yij)t ifte [rj,r,_i] 

for 1 < j < A;, is well defined. Define V by shifting it at the partition points: 

r V{k)t if <e [0,Tfe_i], 

^^^1 V{k~i)t + V{k)T,_, ifte [Tfc_i,rfc_2], 

i ni)t + Et2nOT,_, ifte [Ti,T]. 

Then e C([0,T]; H'*'). Finally we define 

Mt^Yt + Vt fori€[0,T]. 
It remains to show that M is a martingale. 
Lemma 4.2 A/ defined above has the expression: 

k 

= M(j\ + V{1)t,_, if t e [T„ r,_i] (4.7) 
i=j+i 

for 1 < J < fc, and moreover, M is an (Tt) -martingale up to time T , so that 

Mt = E[^ + VT\Tt). 
Proof. We first prove the expression (|4.7p . Since for 1 < j < fc, 
Y{j)t^M{j\^V{])t if [T,,T,_i] 

so that 

k 

Yt - M{])t + V{1)t,_, -Vt if t € r,_i], 
one may conclude that 

k 

Mt = M{j)t + J2 ^(Ot,_, if t e [T„T,_i]. 
i=j+i 

It is clear that M is adapted to (J^t), so we only need to show E{Mt\Ts) = Mg 
for any < s <t <T. If s, t € [T^, Tj-^i] for some j, then 

Mt - M, = MU)t - M{j)s 

so that 

E{Mt - M,\Fs) - E{M{j)t - M{j),\Fs) - . 
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li s e [Ti,Ti^i] and t G [T,, for some i> j, then according to (|4.7p . 

and 

k 

Mt = Mij)t+ J2 nOT,_i. 
i=j+i 

Since M(j) is a martingale on [Tj,Tj^i] so that 

k 

i=j+i 

by conditional on -FTj^i C J^t, we obtain 

fc 

£;(Afi|^T,+J = -E(Af(j)T, + V{j + 1)t,|-Ft,+J + ^(Ot,_,. (4.8) 

On the other hand, M{j)T^ = Yt^ + V{j)Tj = so that 

E{M{j)t, + V{j + 1)t, I J-t,+, ) = E{Yt^ + V{j + 1)t, \Tt,^, ) 

= M(j + 1)t,+,. 

Substituting it into (|4.8p we obtain 

fc 

£;(Mt|J-T,+J = A.f(j + 1)t,+, + 5] ^(/)t,_,. (4.9) 

l=j+2 

By repeating the same argument we may establish 

fc 

E{Mt\TT,_J = M{i - 1)t^_, + Y,yil)T:-i- (4-10) 
Since s G [T^jTi^i], by conditional on J^s, 

fe 

£;(Mt|J-,) - £;(Af(* - 1)t._, + V^(i)T._J^.) + Y ^(Ot,_, 

/=i+l 

fc 

= E{Yt^_,+V{i)t^^AJ's)+ J2 ^Wt,_, 

fc 

= M{i),+ Y ^(Ot,_, 
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which proves M is an J"i-adapted martingale up to T. ■ 

Since L satisfies the local-in-time property and the differential property, so 
that 

i[T,.T,_i](M(VS-))s = L{M)s for s e 

hence 

V{j)t= f h{s,Y,,L{M),)ds + Y. I Ms^Ys)dBl 
for any t e [Tj,Tj^i] and j — 2, - ■ ■ ,k. Therefore 

pt d nt 

Vt^ fo{s,Y,,LiM),)ds + y2 Ms,Ys)dBl Vte[0,T] 
Jo Jo 

and Y = M — V, Yt = which together imply that 

Mt-Yt^ fo{s,Ys,L{M)s)ds + y2 Ms,Ys)dBl Vt e [0,T]. 
Jo Jo 

Thus {Y,M) solves the backward equation (j2.3p . Uniqueness follows from the 
fact the solution {Y{j), M{j) — AI{j)Tj) is unique for any j. 
The proof of Theorem 14.11 is complete. 

We end this article by several comments about the main results. 

The local and global existence results remain valid even if the driver /g 
and the diffusion coefficients of the BSDE are random as long as the global 
Lipschitz conditions are maintained. For example, if /q : i?+ xflx x i?™ 

and ff : i2+ x fi x are jointly measurable such that for any special 

semimartingale F and Z e n^{[0,T];R^) (resp. C{[0,T]; R^)), fi>{t,-,Yt,Zt) 
and f^{t, - jYt) are progressively measurable and 



^ E \ f^it,;Yt,Zt) - f^it,;YuZt)\^dt 

< CallF - i^||c[T2,Ti] + C^WZ - Z||-H2[T2,Ti] 



and 



E r \fl(t, ; Yt) - n(t, -M'dt <Cs\\Y- y||c[T2,T,] 

J To 



for any [T2,Ti] C [0,r], F, f e S{[0,T];R^) and Z, Z e ^^([O, T]; H") 
(and similarly for the case Z,Z e C{[0,T];R^) with norm C[T2,Ti] instead 
of H'^[T2,Ti]), then all our local and global results remain true. We leave the 
details of the proofs for the reader who may be interested in such a generaliza- 
tion. 
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